Low-frequency Noise in Josephson Junctions for Superconducting Qubits 
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We have studied low- frequency resistance fluctuations in shadow-evaporated Al/A10a;/Al tunnel 
junctions. Between 300 K and 5 K the spectral density follows a l//-law. Below 5 K, individual 
defects distort the l//-shape of the spectrum. The spectral density decreases linearly with temper- 
ature between 150 K and 1 K and saturates below 0.8 K. At 4.2 K, the spectral density is about 
two orders of magnitude lower than expected from a recent survey [D. J. Van Harlingen et ai, Phys. 
Rev. B 70, 064510 (2004)]. Due to the saturation below 0.8 K the estimated qubit dephasing times 
at 100 mK are only about two times longer than calculated by Van Harlingen et al. 
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Superconducting qubits arc promising candidates for 
a solid-state realization of a quantum computer. Sev- 
eral groups have demonstrated high-quality operation of 
single or coupled qubits working in the charge Q , charge- 
phase 3: fl^^^ H; or phase regime jj|. Decoherence 
due to external sources such as the measurement devices, 
has been studied extensively and is by now well under- 
stood [3 , permitting qubit dephasing times of up to sev- 
eral microseconds .&j. Future progress in this field of 
research depends crucially on understanding and control- 
ling decoherence due to defects in the devices [3,I1,I3- S^" 
perconducting qubits contain Josephson junctions, whose 
Josephson energy Ej = <I>o/c/(27r) determines the po- 
tential landscape of the qubit {Ic is the critical current, 
and $0 = h/2e is the superconducting flux quantum). 
Due to imperfections of the tunnel barrier, Ej fluctuates 
in time, leading to fluctuations in the qubit potential. 
Therefore, the qubit energy splitting is not constant dur- 
ing an experiment, which leads to decoherence. 

The role of critical current noise in superconducting 
qubits was addressed by Van Harlingen et al. 9]. They 
calculated the sensitivity to /c-noise for different qubit 
designs and estimated decoherence for various data ac- 
quisition schemes, assuming a l//-dependence of the 
noise spectral density. For a quantitative comparison, 
they compiled l//-noise strengths from published data, 
which turned out to be remarkably universal at 4.2 K. 
To extrapolate to millikelvin temperatures where qubits 
are operated, they assumed a T^-scaling of the spectral 
density, based on an experimental study [13 • This as- 
sumption is supported by a recent theoretical paper ^J . 
The experimental data compiled in Ref. |9| did not con- 
tain measurements from the Al/A10a;/Al material sys- 
tem, which is utilized by most groups working on super- 
conducting qubits. Also, no sub-micron junctions were 
included in the survey. For those junctions, it was demon- 
strated that individual bistable defects in the barrier can 
distort the l//-shape of the noise spectral density [l3| . 
Thus, for a reliable estimate of qubit dephasing times, es- 
sential parameters were still lacking. This motivated us 
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FIG. 1: (a) Simplified schematics of the measurement setup 
using a resistance bridge circuit and cross correlation of the 
amplifier signals. (b) Scanning electron micrograph of a 
shadow evaporated tunnel junction, viewing angle 60° to the 
surface normal. The Al leads are widened immediately after 
the junction to enhance cooling. 



to investigate low-frequency /^-fluctuations in Al junc- 
tions prepared by the double-angle shadow evaporation 
method [l3J, which is employed for many qubit experi- 
ments. 

Low-frequency noise of a Josephson junction is caused 
by bistable defects situated in the tunnel barrier |13 |. 
The fluctuating tunnel transparency results in fluctua- 
tions in the tunnel resistance Rn and the critical current 
Ic ■ Each individual defect produces a random telegraph 
signal in the device resistance with a Lorentzian spec- 
trum [illil 
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where SR is the resistance change caused by the defect 
and Toff is the average transition rate between the two 
states. If the random telegraph signals of many weak 
defects with a range of SR and TeS are superimposed, 
the total noise spectrum of the tunnel junction shows a 
1// shape [Tjl- In sub-micron junctions however, a small 
number of strong defects dominate the noise spectral den- 
sity, giving rise to one or more Lorentzians on top of a 
1// background [l2 |. 

Assuming that the /ci?Ar-product is constant, we mea- 
sure fluctuations of the tunnel resistance Rn- This is 
in contrast to most previous measurements, where the 
Ic noise was measured directly, but was shown to give 
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FIG. 2: Noise spectral density of a sample with 880 x 190 nm^ 
junction area from room temperature (topmost trace) to 0.3 K 
(lowest trace) . The data is normalized to the tunnel resistance 
of the sample. The background was determined by measuring 
two metal film 1 kf7 resistors in bridge configuration. 
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FIG. 3; Noise spectral density of the same sample as in Fig. 1, 
for T between 5 K and 0.3 K. At the lowest temperature, a 
Lorentzian spectrum was fitted to the data. 



the same relative noise [Ijl and also allows us to mea- 
sure noise above Tq- Our setup consists of two nom- 
inally identical samples in a bridge configuration [I3 
mounted in a ^He-cryostat with a base temperature of 
about 250 mK (see Fig. ^(a)). The bias resistors have 
more than 1000 times the sample resistance and are sit- 
uated at room temperature. The samples are kept in 
the normal state at all temperatures by applying a mag- 
netic field well beyond 100 mT and are biased with an 
ac current at 1 kHz, which is the noise optimum of our 
amplifiers. We use a pair of home-made battery-powered 
pre-amplifiers at room temperature with an input volt- 
age noise of less than 900 pV/v^Hz, followed by lock-in 
amplifiers. Their output signal is digitized and the cross- 
correlation spectrum is computed to suppress amplifier 
and cable noise. The noise floor of the system is deter- 
mined by current noise injected from the amplifier inputs 
to the sample (about 60 fA/VHz), and residual micro- 
phonic pickup in the voltage sensing cables. The long- 
term stability of the setup permitted measurements down 
to 1 mHz at all temperatures (see Fig. 12). We determined 
the overall noise floor at a typical source impedance by 
measuring two precision metal film resistors of 1 kQ at 
300 mK. 

The choice of bias currents was dictated by heating. 
For test junctions in a geometry similar to flux qubits, 
the poor thermal conductivity of the narrow leads limited 
us to 400 nA. Therefore, we also fabricated samples with 
comparable junction areas, but rapidly widening leads 
(see Fig.n](b)). This permitted bias currents up to 2 ^A. 
Since we suppressed superconductivity in the Al lines, we 
could estimate the thermal conductivity seen from the 
junction from the geometry and the measured sheet re- 
sistance using the Wiedemann- Franz law [I3|- From this, 
we calculated the self-heating to be below 100 mK at the 
highest bias current and at the lowest temperature. We 



also repeated the measurements at different excitation 
levels, to confirm that the noise results were not affected 
by heating. 

We prepared the samples using the double-angle 
shadow evaporation method |l3- The qubit test junc- 
tions were fabricated in the same geometry and on the 
same chip as qubit devices with junctions sizes between 
140 X 250 nm^ and 200 x 250 nm^. The samples with wide 
leads had areas between 490 x 190 nm^ and 880 x 190 nm^ 
and were fabricated in a separate run using the same 
evaporator. 

Figure 121 shows the resistance noise spectral density of 
a pair of two tunnel junctions with 174 fi and 185 fi re- 
sistance and 0.17 /im^ junction area each. The results 
are normalized to the sample resistance, which facilitates 
comparison between different samples. At room temper- 
ature, the spectrum follows a 1// law over the whole 
range from 1 mHz to 60 Hz. The spectral density de- 
creases rapidly when cooling down. Below 5 K the spec- 
trum deviates from the 1// behavior. This is seen more 
clearly in Fig. |31 where we plot the normalized spectral 
density Sr/R^ at low temperature. While the l//-like 
contribution keeps decreasing with T at 1 mHz, Sr/R^ 
saturates around 10 mHz, where a Lorentzian shape dom- 
inates the spectrum. Fitting Eq. ^to the data in Fig.O 
gave Teff = 16 s and SR = 0.85 mil. 

We observed a similar behavior in all the samples 
we measured. At the lowest temperatures, the spec- 
trum consisted of a superposition of a l//-like back- 
ground and one or two Lorentzians in the observed fre- 
quency range. To compare different samples, we plot 
Sr/R^ at 15 mHz 48] of five junctions in Fig. 21 nor- 
malized to a junction area of 1 /im^ assuming a scaling 
of Sr/R^ cx 1/area Q- The data points for the differ- 
ent samples are very close to each other for all temper- 
atures. In fact, the variation between different samples 
is about as large as between different cooldowns (see be- 
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FIG. 4: Temperature dependence of the noise spectral density 
of five different samples at / = 15 mifz. Tlie data was nor- 
malized to a junction area of 1 fim^. Open (o) and solid (•) 
circles are data from the same sample taken in two different 
cooldowns. The star (*) is the average value from Ref. with 
the T'^-dependence (gray line) assumed in that article. The 
solid black line is a fit to a linear T-dependence. 



low). We also plotted the noise prediction from Ref. 9. 
It exceeds our data by more than two orders of mag- 
nitude at 4.2 K [13 • Also, the observed temperature 
dependence deviates considerably from the proposed T^- 
behavior. From room temperature to 150 K, Sr/ E? de- 
creases by about a factor of 100, and then follows a lin- 
ear T-dependence down to about 1 K. From the data 
points for individual samples at lower temperatures we 
find that the noise spectral density always saturates be- 
low a temperature around 0.8 K. The saturation was not 
due to a lack of sensitivity at lower T (cf. Figs.j^lEl) nor 
to heating, but rather due to the saturation of individ- 
ual defects dominating the spectrum, as exemplified in 
Fig. El If many of those defects are present, their spectra 
are superimposed and lead to 1 //-noise with a linear [lj| 
or quadratic ^lOj, ,11] temperature scaling, depending on 
the details of the defect dynamics. This picture breaks 
down if a small number of fluctuators dominate the spec- 
trum. At low temperatures, we assume that their dy- 
namics changes from the thermal to the quantum regime 
and the temperature dependence saturates. In this pic- 
ture, we can assume that the noise level will not decrease 
any further when going to below 300 mK. At 100 mK, 
where qubit dephasing times were calculated in Ref. ^, 
the observed noise spectral density is again comparable 
to their estimate, albeit due to a completely different 
temperature dependence. 

We also evaluated the time constants and resistance 
change caused by the two-level fluctuators. The time 
constants, between 1 s and 60 s, were determined by the 
frequency window of our measurements. Assuming that 
the current distribution in the junctions is homogeneous 
and that a defect can block current flow completely in 
a small part of the junction, we obtain effective defect 
areas of typically 1 nm^ to 2 nm^ . 



The two-level fluctuators were changing with thermal 
cycling. For example, the open and solid circles in Fig.^ 
are measurements of the same sample in two consec- 
utive cooldowns. We observed a Lorentzian superim- 
posed on a 1// background in the first cooldown (o), 
and only the 1// spectrum in the second cooldown (•). 
This leads us to assume that the lowest data points in 
our graph are actually the closest estimate for 1// noise, 
yielding 1.3 x lO^^^/Hz at 15 mHz. This works out to 
2.0 X 10~^^/Hz at 1 Hz, if we assume a true 1// depen- 
dence and normalize to 1 /im^. Note that at 1 Hz and 
at the lowest temperatures the l//-noise approaches the 
sensitivity of our setup. 

For qubit decoherence, we have to distinguish between 
l//-noise and individual bistable fluctuators. In the for- 
mer case, due to the divergence as / ^ 0, the total time 
for an experiment as well as the details of signal aver- 
aging are important, as pointed out in Ref. 9. Using 
spin-echo, the coherence time can be improved only by a 
logarithmic factor Q of the order five. Random telegraph 
noise in the barrier resistance, on the other hand, leads 
to a splitting of the qubit spectroscopy lines and beats 
in the coherent oscillations, if the qubit energy landscape 
switches frequently during an experimental run. Unlike 
l//-noise, it falls off quickly at high frequencies and can 
therefore be efficiently suppressed with spin-echo pulses 
or more advanced pulse sequences, such as series of short 
TT-pulses 123. Dephasing times due to pure l//-noise 
were estimated in Ref. P at T = 100 mK, which gave 
a normalized spectral density of a 1 /xm^ junction of 
8.2 X 10^^'*/Hz at 1 Hz. In our junctions, we observed a 
l//-noise of about a quarter of that value at 300 mK, and 
we assume that the noise has already saturated at that 
temperature. If we suppose that the /c-noise in qubit 
junctions is only due to i?-noise we get twice the dephas- 
ing times estimated in Ref.^. For the three-junction flux 
qubit, charge-phase and phase qubits, this works out to 
1.6 /iS, 3.6 /xs, and 28 /iS, respectively. As an example, 
the spin-echo time of 4 /is in a recent flux qubit experi- 
ment |g] was still limited by the measurement circuitry, 
but already approached the expected spin-echo time due 
to l//-noise (« 5 x 1.6 /is). 

To summarize, we have measured low-frequency resis- 
tance fluctuations in aluminum based Josephson junc- 
tions, as used for superconducting qubits. The noise 
spectral density at 4.2 K is two orders of magnitude lower 
than expected from the literature survey in Ref. 0, and 
we find a linear T-dependence between 150 K and 1 K 
instead of the proposed T^-law. The spectral density 
saturates below 0.8 K, which is due to individual strong 
two- level fluctuators. The dephasing times due to pure 
l//-noise are estimated to be about twice as long as in 
Ref.g 
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